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Abstract 

Ad investigation of three finite- difference methods and their responses 
to the insertion of simulated satellite data is presented, A simple one— level 
' barotropic model is used as the ’’forecast”^ model, while die Mintz— Arakawa 
two-layer model is used to furnish the initial field, the verification fields, 
and the simulated satellite data. The schemes tested are the Shuman, the 
^Matsuno-TASXJ, and an implicit scheme devised by- McPherson. Results 
indicate that the schemes react to inserted data as they would react to un- 
filtered initial fields. The schemes which contain significant implicit viscosity 
are capable of damping the high frequency oscillations which occur after 
insertions} biif. Riich schemes may cause a loss of information. The schemes 
"which contain less damping capability produce ’shock’’ waves which damage 
the forecasts. It is also found that insertion of winds along with temperature 
data improves the forecast considerably. 
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L Introduction 

use of asynoptic satellite data to update integrations of numerical 
r.:atmospheric circulation models has been the subject of much numerical experi- 
"^•Tiientation recently. Most experiments deal with what has come to be known as 
four-dimensional data assimilation, which consists of the insertion of data directly 
into the model at the appropriate grid point and time step of the calculation. The 
'""data are usually atmospheric temperature profiles derived from infrared radiance 
measurements. It is hoped that by replacing enough computed temperature values 
with actual temperature measurements, the other field variables will adjust 
themselves to the observed temperature profiles, and the forecast values will 
approach reality. Simulation experiments have been performed using computer 
produced temperature fields. One major obstacle to the direct insertion of data, 
discovered during these experiments as well as during experiments involving 
actual satellite soundings, is the production by the model of small, "shock" waves 
after the insertion, because of the sudden gradient imposed by the insertion. 

■' These shock waves may eventually dissipate, but they may leave in their wake 
a forecast which hardly shows the effect of the insertion. When actual data are 
inserted, as opposed to model-generated data, the shock effect will probably be 
more noticeable because the model is more attuned to data which it has produced, 
adjustmehrdf the other 1Reld variables is accomplish^ more smoothly. 

In a recent experiment, Jastrow and Halem (1973) have' discovered this to be 


the case. 
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Various means to control these shocks were discussed in the Symposium on 
•Four- Dimensional Analysis at Princeton (1971). Some workers believed that a 
weighting of data before insertion would reduce the shocks while still accomplish- 
ing the same goal of updating the forecast field. Others argued that insertions of 
data with or without weights would accomplish the same results. (Kasahara 
(1972) gives a brief review of most papers presented at the symposium). Re- 
.^-;cently, Hayden- (1973) has shown that by means of "dynamic balancing," the 
-effect of insertions can be enhanced. Dynamic balancing was conceived in an 
effort to deal with the problems of initializing the field of variables for numerical 
models. It consists of adjusting the wind, temperature, and pressure fields so 
that they are in dynamic balance and do not give rise to false waves which may 
disrupt the forecast. The insertion of data is similar in many ways to the initial- 
ization of the field and can thus be treated in a similar fashion. Thus, a localized 
dynamic balancing of the area about an insertion should improve the updating, 
process. 

Experiments and discussions dealing with four-dimensional analysis usually 
contain little mention -of the models and -their~re spec tive differencing schemes. 
Although many experiments have been performed at different research centers 
usii^ different models, there is a need to investigate the effect of the models 
themselves, especially with respect to their characteristic differencing schemes, 
-on conclusions concerning insertions. This study is an attempt to demonstrate 
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the response of various finite-difference schemes to the insertion of unweighted 
data. Using a simple barotropic model as a forecast model and the Mintz-Arakawa 
(M-A) two-level general circulation model as a source of data, three schemes 
will be tested in regard to their ability to assimilate inserted data effectively. 

These schemes are: 1. A Shuman scheme employing centered time differences 
and 9-point spatial smoothing; 2. A Matsuno-TASU scheme, which uses a predictor- 
corrector time differencing and an alternating, staggered space differencing; and 
3, a semi- implicit scheme similar to one proposed by McPherson (1971) using 
much larger time steps than allowable with the first two schemes, 

H. Design of the experiments 

The model used to test the response of the differencing schemes was a 
one-level barotropic model containing three equations relating velocity and 
geopotential height. Data selected for initial conditions, boundary conditions, 
insertion data, and verification data were obtained from the two-layer M-A model 
described by Gates et al. (1971). All the data were from a computer simulated 
February situation. Forecasts were made up to 24 hours at the 500 mb level and 
verified against the M-A fields at 6-hour intervals. The domain for all forecasts 
was limited to an area slightly larger than the North Pacific Ocean, extending 
from approximately 2N to 70N and running westward from HOW to 130E. The 
restricted domain necessitated some specification along the four boundaries. This 
was accomplished by inserting M-A data at the boundary points at each time step. 
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■fliereby preventing boundary effects from ultimately ruining the forecasts. The 

grid system employed was entirely similar to the M-A staggered grid system in 

order to avoid interpolation procedures, although this unrealistically tuned the 

forecasting model to the inserted data. The grid increment in the East- West direc 

tion, AX, is 5 deg longitude, while the North-South uinit increment, Ay» is 4 deg 
1 

latitude. The grid is sta^ered in the sense that not all variables are calculated 
at the same grid points. In the M-A model, the thermodynamic variables are 
calculated on grid points seperate from the points at which wind velocities or 
mass transport terms are calculated. The barotropic model has only three 
variables, but its grid is also staggered, with geopotential heights calculated on 
one grid and the components of wind velocity calculated on a grid staggered 
2-1/2 deg west and 2 deg south of the height grid. 

The equations for the barotropic model are similar to the equations 
governing the free surface conditions of a shallow fluid. The equations contain 
no sources or sinks of energy except at the botindaries, but divergence of the 
velocity fields does exist, and hence vorticity is not conserved. The equations 
are: 


The time increment, ^ t, was six minutes for both the M-A model and 
-4he barotropic model when both the Shuman and Matsuno-TASU scheme 
were applied the semi-implicit scheme allowed for a 30- minute time 
increment. 
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where u is the westerly wind component, v is the southerly wind component, 0 
is the geopotential, f is the corioIis parameter, and m is the map factor. Ihe 
independent variables are the East- West distanccj x, the North- South distance, y, 
and the time, t, 

- The experiments used the barotropic model as a forecast model, while the 
M-A model was regarded as the ideal state or ’’nature’'. By inserting data from 
the M-A model into various numerical approximations to (1) and measuring the degree 
to which the forecast was improved by the insertion, some conclusions could be 
formed regarding the effect of particular numerical schemes on the updating 
process. These insertions were performed at each time step, along paths repre- 
senting the swaths of satellite observations, by merely replacing the calculated 
values found at all grid points within the swath with data from the M-A model. 

The entire domain was covered by the equivalent of three staggered satellite 
passes in approximately 4 to 5 hours of -simulated time and all grid points had 
^4heix -values replaced once and only once in that period. The width of each of 
the three swaths is approximately 800 km, implying a side scan capability of 
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400 km from the sub-satellite point. The inserted data were assumed ideal in 
that no "error*' was added to the M-A data, nor were any grid points "unobservable." 

-For each numerical scheme tested, the experimental procedure was as 
follows: 1. A 24-hour control forecast was made using the prescribed initial 
field. The RMS differences between the forecast and the M-A fields were calcu- 
lated at 4, 5, 6, 12, 18, and 24 hours. The spectrum of kinetic energy over the 
central latitude (39N) was determined at each hour and averaged over 12 hours, 

.>2* '^Using the same initial conditions, a second 24-hour integration was made with 
M-A data inserted at each time step until all grid points were covered once. This 
means that insertions were only performed for the first 4 to 5 hours, while the 
model continued to forecast until 24 hours. The data consisted of all three varia- 
bles (i.e. , u, V, and 0) except in one instance which will be explained below. The 
RMS differences and energy spectrum were calculated for the insertion case as 
well. 3, Another forecast was made differing from the control only in that the 
initial conditions were perturbed by random errors. These errors were uni- 
formly distributed with mean 0. The RMS differences between the perturbed 
-initial field and the non-perturbed initial field were 3. 19 m sec”^ for u, 

1, 02 m sec ^ for v, and 32, 46 gpm for 0g ^ , where g is the gravitational 

-2 

constant, 9, 80 m sec . This forecast was made without any insertions. 

4. Starting \v4ih the perturbed initial conditions, another 24-hour nm was made 
with insertions to update the integration. 
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The RMS errors did not asyxnptote by 24 hours, nor were they expected to 
do so. Yet integrations were not carried out for longer periods because the 
purpose here was to investigate the effect of insertions on short term forecasts 
rather than possible initialization through long term insertions of both synoptic 
and asynoptic data. In f^t, the 24 hour forecast can be regarded as a 19- to . 
20-hour forecast based on an initial field derived after 4 to 5 hours of insertion. 

This 4-step procedure was followed for each of the three schemes tested. 

In steps 2 and 4, however, varying methods of insertion were tested. The details 
and results of these experiments will be discussed in the ensuing sections. 

nL Three schemes 

The Shuman Scheme 

The first scheme investigated is one which is in operation in various 
models, notably at the National Meteorological Center in Washington, D, C, 
Details of this scheme, named after F. Shuman, are outlined by Haltiner 
(1971). Briefly, the scheme uses the eight points contiguous to a central grid 
point in taldng differences (see Appendix). This results in significant spatial 
smoothing at each time step. The time differencing is generally centered except 
that some time smoothing is necessary to avoid the growth of instabilities. 

Implicit viscosity of the Shuman scheme is quite small compared with jmany other 
schemes, so that small waves can grow considerably (see, for example. 
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Grammeltvedt (1969)). It is, therefore, advisable to adjust the initial wind field 
by the use of the balance equation or by eliminating the non-divergent component 
of the winds and substituting the divergent component of the forecast wind field. 

Similar sensitivity was encountered when data were inserted during the ex- 
periments, In addition to small oscillations arising due to physical waves, com- 
putational difficulties were encountered because of the centered time differencing. 
The insertion of raw data at any grid point initiated a computational oscillation 
due to the alternating effect of the centered time differencing. This latter defect 
was remedied by "re-initializing*’ the point of insertion, that is by employing a 
forwai'd time step at the point of insertion and at the time of insertion, similar 
to the forward difference employed over the entire field at the initial time step. 
This simple procedure eliminated the destructive computational mode and reduced 
the error in prediction considerably. Fig; la shows the change in time of u at 
the central grid point for the control forecast, the M-A data, the forecast with 
insertions without adjusting the centered differencing (CEN), and the forecast 
with insertions but with forward differencing applied at the point of insertion 
(FWT>). The insertion takes place for this particular grid point at time step 22 
and its effect on the velocity values is depicted until time step 72, long enough 
for a definite wave pattern to appear. No M-A data, however, were available 
beyond 6 hours (time step 60). The CEN and FWD forecasts differ significantly 
afterwards. The CEN curve indicates a high frequency oscillating wave, while 
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the FWD forecast shows no such tendency. The wavelength of the FWD forecast curve 
seems to be equal to the wavelength of the control forecast curve, although it is 
notably out of phase with the control and is closer to the M-A data. The effect 
of the forward differencing is to have the model "forget" the prior calculated data 
5^T>d begin the computation from an almost new initial state. This effect is even 
more noteworthy when there are random errors present in the initial data. In 
this case, the forward differencing helps cancel not only the effect of the diver- 
gence of the forecast from the M-A model but also the effect of the initial random 
error, as can easily be observed in Fig, lb. The CEN forecast, on the other 
hand, oscillates with even greater amplitude than in the former case because of 
these two effects. 

The growth of RMS errors when initial conditions were unperturbed showed 
the same pattern for all 3 schemes. They all exhibited monotonic increases of 
RMS errors for the first few hours, but with the rate of growth much smaller 
when insertions were performed. Results of the experiments involving perfect 
initial conditions are omitted in favor of the more interesting experiments in- 
volving perturbed initial conditions. Fig, 2 demonstrates the impact of the in- 
sertions on RMS error growth. Starting with initial random errors, there is a 
period of adjustment where there is no significant increase of RMS error even 
for the control forecast, (This adjustment may be related to the filtering of 
high wave numbers present in the initial field which are present neither in the 



11 


forecasts nor in the verification fields from the M-A model. In the real world, 
however, the high wave numbers are present in the verification fields as well 
-and this adjustment period will probably not exist). During this time, the insert 
tions help reduce the RMS errors to levels lower than their initial state. The 
FWD forecast is decidedly more efficient in reducing the error level, so that by 
5 hours, the error level reached is the same as the level of error produced by 
^the FWD forecast starting from perfect initial conditions (not shown). In other 
words, the FWD forecast effectively "forgets” the poor initial conditions and 
bases its predictions on the updated inserted data. 

Matsuno-TASU scheme 

The numerical scheme of the M-A model is one which is temporally a 
predictor-corrector scheme while spatially it uses a time-altematir^ space- 
imcentered fTASU) scheme ^ee Appendix); Arakawa (1972) demonstrates the 
merit of this scheme in a simple example involving the two-dimensional shallow- 
water equations without dissipation terms, but with a point source and point sink 
present in the region. All other schemes tested showed a false alternating 
pattern instead-of a monotonic decrease, of the height toward the sink. The 
>Matsuno-TASU scheme, however, produced results which were realistic with a 
decrease in height from every direction towards the sink. The schemers effec- 
tiveness is due to the combination of the Uvo-step Matsimo scheme and the alter- 
nating space differencing of the TASU scheme. 
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The Matsimo schenae involves predicting an intermediate value of the varia- 
ble being dealt with, then using that value to arrive at the desired forecast. That 
is, if Q is the variable to be predicted at time step n+1, given Q at time n, then 
a quantity Q* is determined so that 

Q* = + f(Q°) , (2) 

where the superscript refers to the time level and f^) is the algorithm involving 
^known values by which the prediction is made. Q* is now substituted for in 
the algorithm to obtain 

= q“ + f(Q*) (3) 

The Matsuno scheme, although comparatively costly in computer time, can 
because of its implicit viscosity, filter some of the small waves efficiently. This 
viscosity, after insertions are made, may diminish their beneficial effect while 
damping the shocks. This aspect will be discussed more fully below. 

The TASU scheme involves shifting the space differences to alternately the 
upper-right then the lower- left portion of the grid box surrounding any point. In 
the e^eriment here and in the M-A model the TASU step is invoked only once 
every five time steps and then only during the second half of the time step. At 
all other times, the space differencing is centered. 
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In the experiments involving the Matsimo-TASU scheme, insertions of two 
kinds were made. Since the Matsuno time steps consists actually of two half- 
steps, i. e, one predictor and one corrector, it is possible to insert values after 
either half of the time step. That is to say, it is possible to substitute at a par- 

ticulax grid point at some given time an observed value of Q for either Q* or Q 

n+1 

. . (where Q, Q*, and Q have the same meanii^ as in (2) and (3)), Upon attempt- 
ing both forms of insertion, it was foimd that inserting after the first half step 
given by (2) (herein referred to as ’’insertion after 1”) was far inferior to insert- 
ing after the second half given by (3) (herein referred to as "insertion after 2"). 
This was probably due to the strong damping effect the second half step had on 
the inserted data, their influence being greatly reduced and the forecast. barely, 
if at all, improved because of the insertion. 

Pig, 3 demonstrates the effect of the insertions on the RIVK error for the 
experiment where initial error was present. Only insertion after 2 results in 
a meaningful reduction of RMS error. Insertion after 1 is almost indistinguish- 
able from the control case. All cases resulted in an initial adjustment period 
during which the RMS error was substantially reduced. When insertions took 
. - place during the first few time steps, the RMS errors were reduced to values 
significantly lower than the RMS errors of the control forec^t. The appearance 
of this adjustment period and its obvious reinforcement by the insertions may 
possibly be explained as the result of filtering by both the barotropic and the 
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M-A models. By perturbing the initial conditions with random error, small 
waves are probably produced. The M-A model successfully filters these waves, 
so that the verification fields and the satellite data are devoid of these high wave 
= numbers. The barotropic model because of its spatial resolution and with the 
assistance of the Matsuno-TASXJ scheme begins filtering these waves as well, 
maldng its early forecasts more compatible with the M-A fields despite the pre- 
dictability factor. If M-A data are inserted, -the filtering is apparently acceler- 
ated and the RMS error reduced finally to a level which is nearly independent of 
the initial state. Whether this will occur if actual atmospheric data are used is 
doubtful, but can only be verified through actual experiment. 

Fig. 4 displays the time history of u at the central point of the grid network. 
It demonstrates clearly the effectiveness of the insertion on updatir^ the model. 

. -The control case and insertion after 1 result in values which are substantially 
different from the M-A data. When insertion after 2 takes place, however, the 
value of u is rapidly brought closer to the M-A data and stays there for some 
time. None of the severe oscillations produced by the Shuman scheme are 
present in this experiment, indicating-once again the strong damping of the 


Matsuno-TASU scheme. 
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Semi- Implicit Scheme 

Implicit and semi-implicit schemes axe often more economical than explicit 
schemes subject to the Courant— Friedrichs— Lewy condition for linear stability. 

By expressing more than one term at the forecast (n+1) time level, larger time 
increments are possible without causing instabilities. The larger time steps, 
however, imply that longer periods can be considered synoptic with respect to 
- ^the “insertion of data. Satellite data which are gathered over the entire period of 
the time increment At, are inserted at one instant. What effect this wholesale 
insertion has on the updating process depends on many factors, among them the 
nature of the data themselves and their change with time. 

The experiment conducted involved a scheme recently developed by 
McPherson (1971). This scheme separates the advective terms from the other 
terms of the equation which are evaluated at the forecast time (see Appendix 
for details). By cross differentiation the velocities can be eliminated, leaving 
a Helmholtz equation of the form: 


where 


^n+1 , ^n+1 ^n -1 , _n ^^2 m-1 
V 0 — \0 (D 0 ) = X F. . (D 0 ) 
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( 4 ) 


2 At, n^l 


" 2 2 
X = 1 + f D 


f 



16 


and F is some known function of space at any given time step. V is a numerical 

2 

approxixnatioii to the Laplacian operator V . The superscript represents the 
time step at which the variable is being evaluated. Once 0 is evaluated from 
^), it is a simple matter to determine u^^^ and v^^^ . is solved by using a 
relaxation technique, starting from some initial guess of 0*^ ^ . For this set of 
experiments, a time increment of 1/2 hour was selected, which reduced compu- 
tation time to about 1/2 of that required for the Matsuno-TASU scheme despite 
the iterations required for the relaxation. 

In addition to the insertion experiment involving the assimilation of data 
from all variables at each time step, an attempt was made here to measure the 
efficacy of inserting geopotential heights alone. This is similar to the practice 
of inserting only temperature l^se rates which are derived from infrared radi- 
ance measurements. When contrasted with results obtained from the insertion 

of all variables, the experiment could indicate whether the "balancing” of the 
velocity fields is substantially beneficial to the assimilation process given 
the grid size and the data characteristics for the experiment. When perfect 
initial conditions are used, the insertion of all variables does not create 
a significant improvement over the insertion of the temperature field alone. There 

is, expectedly, a lower level of error when all variables are inserted simply be- 
cause the quantity of inserted data is greater. When random error is added to 
the initial conditions, however, the dynamic balancing accomplished by the inser- 
tion of wind data along with the heights produces significant improvements in the 
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forecast. Fig. 5 shows the R3VIS error growth for u when random errors are 
present in the initial state. As in the case of the Matsuno-TASU scheme, the 
first few time steps result in an adjustment which lowers the RMS error even 
if no insertions are performed. The degree of adjustment is limited when only 
heights are inserted, but a significant dip occurs when all variables are inserted. 
The level of RMS error reached through insertion at 5 hours is again equivalent 
to the level reached when inserting after perfect initial conditions O^ot shown) 
indicating that the insertions update the computations efficiently despite the 
larger time increment of the semi-implicit scheme. 

Fig. 6 depicts the time change of the variable u at the central grid point 
for the various experiments with the semi-implicit scheme. Because of the 
larger time increment, only 15 time steps are presented, so that the simulated 
time lapse is consistent with earlier figures. Etespite the random errors present 
in the initial conditions, the control run does not give rise to oscillations. When 
alone is inserted, the change in u values is not significant for at least 15 time 
steps. The insertion of all variables does, however, have a notable effect which 
-gives rise to oscillations similar to, but not quite as large as, those produced 
- by the Shuman scheme. The large time steps are perhaps a factor in creating 
these oscillations. Large doses of data are injected at each time step and only 
limited filtering can be performed because of the few steps that follow the 
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insertion. It may prove beneficial to reduce the time step immediately 
foUo’wing insertion so that the adjustment to the model may take place in a 
relatively short period of simulated time. 

IV# Energetics and Comparisons 

A study of the energetics involved in the various numerical schemes was 
- performed by calculatii^ the kinetic energy over the central latitude as a 
function of zonal wave number. The wave numbers mentioned here are not 
global but refer to the domain of integration. Only one third of the Northern 
Hemisphere is represented on this domain, so that wave numbers should be 
multiplied by 3 in order to correspond to full- latitude wave numbers. The 
kinetic energy was calculated at each hour and averaged over 12 hours. Only 
one latitude was involved; too much spatial averaging may have resulted in 
disguising the effect of the insertions. 

When perfect initial conditions were used, few departures from the con- 
-v.trol energy spectrum were observed after insertions. When random errors 
were present initially, however, the insertions managed to return the spectra 
to the non-perturbed control spectra. The effect of the random errors was to 
shift energy from low wave numbers to higher ones probably because of the 
-creation of small wave disturbances. The insertion of data from the M-A 
v:model counter-acted this trend -and shifted the energy back to the imperturbed 
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spectral distribution. Fig. 7 shows the various spectral distributions with and 
without insertions. For all 3 numerical schemes, the optimum method of insert 
tion (from an RMS standpoint) is seen to increase kinetic energy in the region of 
wave numbers 1 and 2. Between wave numbers 2 and 4 there is a crossover, 
and kinetic energy is lower for the insertion experiments for these higher wave 
numbers. Wave No. 5 is anomalous in all three experiments, a local minimum 
beii^ achieved by the control but not by the insertion spectra. Beyond No. 5, 
there are variations among the three schemes, but this is to be expected in the 
sensitive low energy regions. 

The 3 schemes also demonstrated varying degrees of skill in forecasting. 
The Matsuno-TASU scheme reached lower levels of RMS error than either of 
the two other schemes. The Shuman scheme, although greatly enhanced by 
the insertions, has a large error growth'rate. The semi-implicit scheme 
produces RIVIS errors at 24 hours which lie between those produced by the other 
2 schemes. No generalizations should be drawn from this behavior, however, 
as to the forecasting ability of a particular scheme. It must be remembered 
that the data for verification were derived from the M-A model, which itself 
-uses a Matsuno-TASU differencing scheme, and not from observations of the 
real atmosphere. All that can be properly assumed is that data produced by 
the barotropic model using the Shuman scheme do not correlate well with data 
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produced by the M-A model using the Matsuno-TASU scheme, A better test of this 
would be to use the same model, but with different finite-difference schemes, 
to produce forecasts which could then be examined and compared. 

V. Conclusions 

Three numerical schemes have been examined with respect to their reac- 
tion to the sequential insertion of data into various grid points during computa- 
tion. Comparisons were made on the basis of RMS errors , the time history at 
an individual grid point, and the 12-hour spectra of kinetic energy. 

Results indicate that inserting can be an effective means of updating model 
integrations regardless of the numerical scheme, but that some care must be 
taken as to how the insertions are to be performed. The Shuman scheme proved 
very sensitive to the inserted data, and severe oscillations could result if cen- 
tered time differencing is retained during the insertion process. The Matsuno- 
TASU scheme, on the other hand, proved capable of filtering spurious short 
waves effectively, but insertions had to be made after the complete time step 
or the influence of the data might be lost. The semi-implicit scheme gave rise 
to some oscillations despite its damping mechanism, probably due to the larger 
dose of data inserted at each time step and the smaller number of time steps 
needed to achieve a forecast. In general, the schemes behaved with respect 
to insertions as they usually act with respect to initial conditions. The Shuman 
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-scheme requires some adjustment of the initial fields (see for example 
'Haltiner (1974, p, 269) ) because of its lack of internal viscosity. The 
r-Matsimo-TASU scheme on the other hand requires no balancing of the initial 
afield (see Somerville, et al, (1974)), in all likelihood because of its strong 
filtering and damping properties. These schemes behave entirely similarly 
with respect to data insertion, as described above. 

The impact of data insertion on tlie energy spectra seemed to be contrary 
to the expected results. Instead of creating small shocks which give rise to 
high- wavenumber oscillations, the insertion of computer-produced data 
transferred high- wavenumber energy created by poor initial conditions to the 
lower wavenumber rcgidn. This may be seen as due to essentially the re- 
initialization of the entire field with the inserted data, which have already been 
carefully filtered by the M-A model. If real atmospheric observations were 

used, this shift in the energy spectra may well have been in the other direction. 

1 

hi the future, research should center about the use of actual data rather 
than computer simulations. This should help determine which of the numerical 
schemes are best equipped to handle expected satellite observations of the real 
atmosphere and what possible methods can be introduced to improve their 
assimilation. 
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Appendix 


Details of the schemes 

1* Shuman scheme ~ The details of the Shuman scheme are de- 
cribed in Haltiner (1971). (1) written in terms of the scheme become: 
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differences between the equations presented here and those given by Haltiner 
(1971, p. 228) because of the staggered grid system employed in these experi- 

i 

ments. 

All terms are evaluated in the center of a grid box consisting of 4 points. 
The final value for a grid point is determined by averaging the values in the 
four surrounding box centers. Extensive smoothing is thus accomplished at 


each time step. 
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2* Matsuno-TASU scheme - A complete description of this scheme is given by 
Arakawa (1972), The Matsimo time differencing scheme involves a predictor- 
corrector computation with two iterations at each time step as described in 
section 3* The TASU space differencing alternates wdth centered differencing 
once every 5 time steps. The TASU differencing involves taking x-differences 
centered Ay/2 above the central point and ax/2 to the right for y- differences 
during one time step, then taking differences below and to the left for the next 
time difference. The TASU differencing is performed only during the second 
of the Matsuno time steps. The sequence will then look like this for six time 
steps; 


Time step 

Time difference 

Space difference 

1 

Forward 

Centered 


Backward 

Centered 

2 

Forward 

Centered 


\ 

Backward 

Centered 

3 

Forward 

Centered 


Backward 

Centered 

4 

Forward 

Centered 


Backward 


Centered 
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Time step 

5 

6 


Time difference 
Forward 
Backward 
Forward 
Backward 


Space difference 
Centered 
Upper-right 
Centered 
lower-left 


3, Semi- implicit scheme - Many implicit and semi-implicit schemes have 
been devised to avoid the restriction on the time step imposed by the 
Coixrant- Friedrichs- Lewy condition for linear stability. The one used for this 


study is based on the scheme devised by McPherson (1971). 

The equations of motion are written in difference form with all but the 
advective terms given at the same time step; 


- u^)/D +m(^u 
,n+l 


. n ^ n 

n 3u n Bu 

+ V 


bx 


dy 


) 


+ m 




Bx 


- n+1 

- f V = o 


. n+1 k. /-._n 5v . n Sv ^ . 50 . . - 

(v -* V )/D + m ( 'i — + v + f u 

^ V 5y y dy 


n+1 


= o , 


where D = 2 At and k = n-1 except when n = o. (At n = o, D = At and k = o.) 
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Eliminating from the first equation and from the second equation 
leaves explicit equations for and in terms of as follows: 


, n+1 „ ^ 

Xu - B - D f m 


30 


n+1 


n+1 


30 


3x 
n+1 


+ f D 


30 


3y 

n+1 


Xv*‘* = A-Df 

V 3y 


- f Dm 


30 


3x 


) 

), 


where X = 1 + 




A = v -fDu + D 


^ f D u 


(A2) 


k _ r - _ n 3u ^ n 3u 

3x 3y 


n 

r) 


- u 


n 3v n dv 

+ V 


n 


^ J , and B=u^ + fDv^ ~D fD 


dx 


dy 


^ n n n 

n 3v X n Su n 3u 'n -i 

1 ^ ^ ✓" N« ^ J ^ _i 


Cross differentiating and adding to find the divergence which can then be 
substituted into the geopotential equation to give; 


„2 ,n+l , ^n+1 _2 ^n ,~1 , _n _2 n -1 

70 -X0 (D0) =XFp0) 


^3) 


2 

2 

3 

3 

2 

2 

3X 

3y 


V 3y 


3x ay m —-0 IS a known 


3y 


function of space at each time step. 
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(A3) is the equation which appears in section 3 and can be solved numerically 
by means of a relaxation technique. The relaxation was performed sequentially so 
that lower-left values were immediately replaced by the new values. The remainder 
: R was thus determined at each point jj't by the formula; 


^m+1 2 , 2 , -2/'^m , ,n+l > 


(A4) 


- (2 + 2 ( Ay)^ ( Ax) ^ + X ^ ) ~ X 0^ ) ^ , 


where the superscript in terms of m refers to the iteration number. Once R 

n+1 

is determined, the value of 0 is calculated by over-relaxation, i.e, adding to 

.m.._ _ r ^ ^ 2 _ ,2 _ -2 , ^2^n,-li 

0 the quantity i, ok + ( Ly) { Ax; -t a. \u } j • 

-3 

continued until all R , ^ lO , The number of iterations required for this ac- 

j.-t- 

curacy ranged between 35 and 40 at each time step. 
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Fig. 1: 


Fig, 2: 




Fig. 4: 


Fig. 5: 


Figures 

a) Time change of u (m sec at the central grid point using the 
Schuman scheme, for the control run, the insertion run maintaining 
centered differences (CEN), the insertion run with forward differences 
at the time and point of insertion (FWD), and the M-*A data to 6 hours. 
The time step is 6 minutes, b) Same as a) except that random errors 

are present in the initial conditions. 

-1 

Growth of error in u (m sec ) for the Shuman scheme for the 
control case, insertions with CEN, and insertions with FWD. Random 

errors were present in the initial conditions. 

GrovAh of error in ii (m sec for the ^ats\,ino“TASTJ scheme 
for the control case, insertions after the first half time step 
(after 1), and insertions after the full time step (after 2). Random 

errors were present in the initial conditions, 

-1 

Time history of u (m sec ) at the central grid point during Matsuno- 
TASU experiments for the control case, insertion after 1, insertion 
after 2, and the M-A data to six hours. Random errors were present 
in the initial conditions. The time step is six minutes. 

Growth of Ri\iS error in u (m sec“^ ) for semi- implicit scheme for 
the control case, insertion of all variables (u, v, and 0), and inser- 
tion of 0 alone. Random errors were present in the initial conditions. 
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Fig. 6; 


Fig. 7: 


Time history of u at the central grid point during semi- implicit 
scheme experiments including the control, insertion of all varia- 
bles, insertion of 0 only, and the M-A data to six hours. Random errors 

were present in the initial conditions. The time step is 30 minutes, 

2 -2 

Kinetic energy (m sec ) spectra averaged over 12 hours for the 
central latitude ^9N) for the control case and the most effective 
insertion of a) the Shuman scheme b) the Matsuno-TASU scheme 
and c) the semi-implicit scheme. 
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